AN ARTIFICIAL VISCOSITY APPROACH 
TO QUASISTATIC CRACK GROWTH 
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Abstract. We introduce a new model of irreversible quasistatic crack growth in 
which the evolution of cracks is the limit of a suitably modified e-gradient flow of 
the energy functional, as the "viscosity" parameter e tends to zero. 
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1. Introduction 

In this paper we consider the crack growth in brittle materials in the particular case 
of a preassigned crack path Y. We assume that T is a regular arc with one endpoint on 
the boundary of the reference configuration Q, that there exists an initial connected 
crack starting from the boundary point, and that the crack remains connected during 
the evolution. Hence, such a crack will be completely determined by its length a. 
The evolution is assumed to be irreversible, so that the length of the crack will be 
increasing in time, and quasistatic, i.e. at each time the configuration describing the 
body is in equilibrium. By configuration we mean a pair (u, a) where u represents 
the displacement, and a the length of the crack. 

The main new feature of this model is that it is based on a local stability criterion 
for the energy functional rather than on a global one. The choice of the total energy 
£(t)(u, a), of a configuration (u, a) at time t, is inspired by Griffith's idea [11] that 
the evolution of cracks in brittle materials is the result of the competition between the 
elastic energy of the body and the energy needed to extend the crack. For simplicity 
of exposition, we detail our study in the case of antiplane shears and therefore, u is 
a scalar function representing the displacement orthogonal to the plane of Q and the 
bulk part of the energy is given by the square of the L 2 -norm of the gradient of u. 
We consider the case of a homogeneous isotropic material and, according to Griffith's 
theory [11], we assume the surface energy to be proportional to the length a of the 
crack, the constant of proportionality being given by the toughness of the material. 
The evolution is driven by time-dependent imposed boundary displacements ip(i) on 
a part 8d& of the boundary, and applied boundary forces g(t) on the remaining part 
<9at£7. The total energy, £(t)(u,o~), is the sum of the bulk energy and the surface 
energy minus the work of the applied forces g{t). 

Let AD(ijj(t), a) denote the set of admissible displacements, i.e. displacements with 
finite bulk energy, compatible with the imposed boundary displacement ip(t) and with 
the crack length a. Let us recall that the functional S(t)(u,a) is not differentiable, 
nor convex. It depends on a both through the surface energy term and through the 
constraint on the set of admissible displacements. Note that, given t and a, there 
exists a unique minimizer u t ^ of the energy £(t)(u, a) in AD(ip(t), cr). Then let us 
consider the minimal energy corresponding to the boundary data ip(t) and to the 
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crack length a: E(t,a) := £ (t)(u ttC7 , cr). The derivative d a E(t,a) can be computed 
(see Proposition 2.4) and it is related to the stress intensity factor of the displacement 
u tt a at the tip of the crack. It plays a crucial role in the Griffith's criterion for the 
propagation of cracks. 

Let us define now the notion of evolution we are interested in. The irreversible 
quasistatic evolution problem consists in finding a left-continuous function of time 1 1— > 
(u(t), cr(t)) such that the displacement u(t) at time t belongs to the set AD(ip(t), a(t)) : 
and the following three conditions are satisfied: 

(a) local unilateral stability: at every time t > 

£(t)(u(t), a(t)) < £(t)(v, a(t)) Wv G AD(ip(t), a(t)) 
d a E(t,a(t)) > 0, 

(b) irreversibility: the map t h- > a(t) is increasing; 

(c) energy inequality: for every < s < t we have 

£(t)(u(t), a(t)) < £(s)(u(s), a(s)) + Work(w; s, t) , 

where Work(w; s, t) denotes the work of external forces. 

A solution to this problem will be called an irreversible quasistatic evolution. 

We show that conditions (a)-(c) are enough to ensure that at almost every time t a 
weak version of Griffith's criterion is satisfied (see Proposition 2.6). In condition (c) 
two terms contribute to the work of the external forces: the first is due to the surface 
forces generated by the imposed boundary displacement and the second one comes 
from the applied surface loads. 

In [9], [4], [2], [8], [3] the globally stable irreversible quasistatic evolution problem 
was studied. It is a particular case of the previous one and it fits the general scheme of 
the continuous-time energetic formulation of rate-independent processes developed by 
Mielke and his collaborators (see [17] and the references therein). It consists in finding 
an irreversible quasistatic evolution which satisfies the global stability condition: at 
every time t > 

£{t){u{t),a{t)) <£{t){v,a) V a > a(t) V v G AD(ijj(t), a) . 

In this case condition (c) can be replaced by the energy balance condition: the in- 
crement in stored energy plus the energy spent in crack increase equals the work of 
external forces. 

The global minimality condition imposes the comparison, in terms of energy, of 
a configuration with all admissible configurations with a longer crack and might 
generate jumps in the length of the crack that are not justified by the mechanical 
interpretation of the problem. That is why we look for a selection criterion different 
from the global stability. 

The selection criterion we choose is based on an approximation procedure with 
a regularizing effect. We consider solutions (u £ (t),a £ (t)) of regular evolution prob- 
lems that converge to a solution (-u(t), <x(t)) of the irreversible quasistatic evolution 
problem. Moreover, among the possible approximations we choose one that has the 
following property: 

(V) if on a certain time interval [^1,^2] there exists a regular function o- (t) such 
that 



d a E(t, a (t)) = and 8 2 a E(t, a (t)) > V* G t 2 ] , 
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and if & £ (t) > for every t G [£1,^2], then the equality o"(ti) = o"o(^i) implies 
that tj(t) = cr (t) for every t G [ti,t 2 ]- 

Let us now describe more in detail the construction of the approximating evolutions. 
First of all, we fix an initial condition: assume that at time t — the crack length is 
equal to a > and the displacement is equal to Uq, in such a way that the initial 
configuration (uq, <j ) is in equilibrium. Then, for every a between o"o and a, where a 
is the length of T, we consider a diffeomorphism $ CT of Q that transforms the crack of 
length a into the one of length a . Using we change variables in the expression 
of the energy functional £ and transform it into a functional T depending on the 
time t, the crack length a, and the modified displacement v, which, assuming that 
the change of variables does not influence the terms due to the boundary data, has 
the form 

F(t,v,a)= / (A(a,x)(Dv(x) + Dip(t,x))\Dv(x) + Dip(t,x))dx + a - 
Jn\r(a ) 

- g(t,x)(v(x)+ip(t,x))dH 1 (x). 

Here, r(<7o) is the crack of length <7o, A(a,x) is a 2x2 symmetric matrix of smooth 
coefficients coming from the change of variables, Dv is the distributional gradient 
of v with respect to the spatial variables x, (-|-) denotes the scalar product in M 2 , 
and Ti, 1 the one-dimensional Hausdorff measure. The advantage of this change of 
variables is that now the set of admissible functions v does not depend on t, nor on 
the crack length a. The same change of variables is considered, in a suitable small 
neighbourhood of the crack tip, in order to compute the derivative d a E(t, a) (see also 
[13], [1], [15]). 

The approximating evolution (v £ ,a £ ) is the solution of a suitably modified e- 
gradient flow for the functional T which starts from the initial data (wo? o"o) : 

(ev £ = -grad^t,^,^) 
e& £ = (-9 CT J r (t,t; £ ,cr e )) + A(o- e ) 
v £ (0) =Uq 
(j £ (0) = (Jo . 

Here grad w JF(t, v, a) denotes the gradient of the function v 1— > jF(t, v,a) considered 
as a function defined on the Sobolev space H 1 (£l\T(a )) with suitable boundary 
conditions. The positive part in the second equation guarantees the irreversibility of 
the evolution, while A is a Lipschitz continuous positive cut-off function that becomes 
zero for a = W, so that only increasing solutions with crack length less than a are 
considered. If we are interested in the evolution until a certain crack length a±, with 
o"o < Ci < is reached, then we choose A(er) = 1 for cr < a < <j\. In this way 
for crack lengths less than <Ji, the regularized evolution law is proportional to the 
gradient flow for JF, with the constraint that the crack length is increasing, while 
it is distorted by A for crack lengths between <j\ and a. Therefore the evolution is 
considered meaningful only until the crack reaches the length <j\. 

Note that, using the form of the functional J 7 , the first equation in (1.1) can be 
written as 



eA x v £ (t,x) = -div x (A(a £ (t),x)Dv £ (t,x)) + . 
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with suitable boundary conditions. We preferred the evolution problem in H to the 
usual parabolic one 

ev £ (t, x) = -dxv x (A(<j £ (t),x)Dv £ (t, x)) + . . . 

which corresponds to the gradient flow in L 2 , because it helped us to prove prop- 
erty (V), see Theorem 4.1. Note also that in this way the first equation in (1.1) 
becomes an ODE and thus the existence of the solution for this modified e-gradient 
flow follows from classical existence and uniqueness results for ordinary differential 
equations in Banach spaces. 

Let us remark that this model is not suited for the study of the crack initiation 
problem. We also note that the approximating evolutions we consider have been 
chosen on the basis of their mathematical simplicity and do not seem to have any 
mechanical interpretation. Nevertheless, we think that the notion of approximable 
irreversible quasistatic evolution proposed here could be the starting point for the 
study of different approximations with a mechanical justification. For a different 
approach to the irreversible quasistatic crack growth see also [10]. 

In Section 5 we detail our results in the case of monotonically increasing in time 
imposed boundary displacements and compare this evolution with the one proposed 
by Francfort and Marigo in [9], while in Section 6 we provide an example where the 
energy, as function of the crack length, has at least a concavity interval. 

2. Setting of the problem 

2.1. The reference configuration and the crack. Let Q, be a bounded connected 
open set of 1R 2 with Lipschitz boundary dfl. The set f2 represents the reference 
configuration of an isotropic, homogeneous elastic body. Let doQ be a closed subset 
of dQ with ^(cfefi) > 0, where Ti 1 denotes the one-dimensional Hausdorff measure, 
and let d^^l '■= dfl\do^- On the Dirichlet part of the boundary, do^, we will impose 
the boundary displacements, while on the Neumann part of the boundary, djyQ, we 
will prescribe the boundary forces. 

Let T be a simple C 3 -arc and let 7: [0,cf] — > T be its arc-length parametrization. 
We assume that 7(0) G c^fi and 7(d) € Q for < a < a. For technical reasons 
it is convenient to extend T until it reaches another point in <9jvf2, so that it cuts 
the reference configuration f2 into two subsets. The extension will still be called T, 
and its arc-length parametrization will now be 7: [0, c max \ — > T. We assume that its 
intersection with the boundary dVt is not tangential. Let v be a unit normal vector 
field on T. Then we denote by Q + the part of \ F which is positively oriented 
with respect to v, and by VL~ the remaining part, so that f2 \ T = f2 + U fi~. Both 
Q + and Q~ are bounded connected sets with Lipschitz boundary. We assume that 
H l {d D Vt n dtt + ) > and H\d D Q. H dQr) > 0. We make the following simplifying 
assumption: all admissible cracks are of the form 

r((j) := {7(s) : < s < a} witha<a. 

According to Griffith's theory we assume that the energy spent to produce the 
crack T(a) is proportional to the length of the crack, and, for simplicity, we take it 
to be equal to a. 

2.2. The bulk energy. We consider here the case of antiplane shears. Given a 
crack T(a), an admissible displacement is any function u E if 1 (Q \ T(a)), and the 
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bulk energy associated to the displacement u is 

W(Du) := [ \Du(x)\ 2 dx, 

where Du is the distributional gradient of u and | • | denotes the norm in M 2 . 

2.3. The boundary displacement. In the following it will be convenient to work 
on a fixed time interval [0, T] with T > 0. We impose a time-dependent Dirichlet 
boundary condition on dr)Q: 

u = ip(t) on do^i, 

where equality on the boundary is considered in the sense of traces. We assume that 
ip(t) is the trace on d D VL of a bounded Sobolev function, still denoted by t i— > ip(t), 
with ip(t) G -ff 1 (o) n L°°(n). 

We assume also that j> G W^°°{0,T; H 1 ^)) n L°°(0, T; L°°(fl)). Thus, the time 
derivative t h- > tp(t) belongs to the space L°°(0, T; and its spatial gradient 

t i-> Dij>(t) belongs to the space L°°(0,T; L 2 (fi; M 2 )). 

2.4. The external loads. We are interested in the case of time-dependent dead 
loads, in which the density, g : [0, T] x d^tt — > R, of the applied surface force per 
unit area in the reference configuration does not depend on the displacement u. We 
assume that the function t h- > •) belongs to W 1,oo (0, T; L 2 (8n^I, 'H 1 ))-, with time 
derivative denoted by t i— > g(t,-). The associated potential, for a displacement u, is 
given by 

0(i)(u) := / g(t,x)u(x)dH 1 . 

Moreover, assume that for every t G [0, T] the support of g(t, •) does not intersect the 
set r. 

2.5. The admissible displacements and their total energy. For every t G [0, T], 
the set AD(ip(t), a) of admissible displacements in Q with finite energy, corresponding 
to the crack T(a) and to the boundary data ip(t) is given by 

AD(iP(t), a) := {u G H\Q\ T(a)) : u = ip(t) on d D n} , 

where the last equality refers to the traces of u and ip(t) on d D Vl. The total energy 
of a configuration (u, a) with u G AD(ip(t), a) is given by 

5(t)(u,<r) :=W(Du)+<t-^(*)(u). 

Note that it does not depend on the particular extension ip(t) chosen, but only on its 
value on the Dirichlet part of the boundary. 

2.6. Moving to a fixed domain. Let Hq n (Cl\T(a)) denote the space of functions 
u G H 1 ^ \ r(er)) whose trace on dpVt is zero. We may consider the energy as 
a functional defined on Hq dQ (Q \ T(a)) by simply writing u = u + ip(t) with u G 
AD(ip(t)i<r) and u G HQ Dn (il \ T(a)). Still the domain of the functional would 
depend on a. To transform it into a functional defined on a fixed domain we consider 
the following change of variables. 
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For a G [o"o,cr], let $(-,cr) = $o-(-) : — > be a diffeomorphism which coincides 
with the identity near the boundary of fi, leaves invariant both f2 + and Q~ and 
transforms r(<r) into T(a ). Let *(-,<r) = := $ _1 (-,(r): ft -> ft. Then 

/ |£>u + £ty(i)|fo; = / | J D M (^(t/)) + J D^(t)(^(y))| 2 det J D^(t/)^. 
^\r(a) Ja\r(a ) 

For « G if^ n (ft\r(<7)) define := u(V a (y)) and let #)(y,<r) := ^(t)(V a (y)). 

With these notations 

/ \Du+D^(t)?dx= [ \((D^ a fr\y)(Dv(y,a)+Dm(y, ( r))\ 2 detD^ (T (y)dy, 
Jn\r(a) Jn\r(a ) 

and the last integral can be written also in the form 

/ a ij(y^ a ) D j( v (y^ a ) + 4>( t )(y,< T )) D i( v (y, (T ) + 4>(t)(y,<?))dy, 

Jn\r(a ) lJe{li2} 

with the coefficients a^- given by the change of variables. 

Define A(a) := (o>ij(cr))ij and note that aij(a) G C(ft), and a,ij(a) = aji(a), for 
every a G [er ,cf], and every 

We may assume that < c < Hdet-D^o-Hoo < C independently of a G [00, ef], where 
II " II oo denotes the L°°-norm on ft. Since T is of class C 3 , we may also choose $(•, <r) 
(and hence to depend regularly on a in such a way that, as functions of a, 

the coefficients be of class C 2 on [<7o,<7], uniformly in ft. In particular, we shall 
use the fact that there exist positive constants A, A, A', L , L' > independent of a, 
such that 

(A(a)^)>\\^\ 2 V£gM 2 , Vxefl, (2.1) 
where (-|-) denotes the scalar product in R 2 , 

IKA^WIU^A^IH V^gM 2 , (2.2) 

IK^^WIU^AleiM V^^GM 2 , (2.3) 

\\a l3 (a') - Oij-Mlloo < - ^"1 and (2.4) 

Rc^V) - ^-(Olloo < - A (2-5) 

for every a', a" G [o" , (f] and i,j = 1, 2. 

Note that, since coincides with the identity near the boundary of ft, this change 
of variables does not have any effect on Q: 

Moreover, we can neglect the dependence of ip on a since, for every a G [oojZf]) ^ a 
coincides with the identity near the boundary of ft, and we may assume that the 
support of ip is included in the set where, for every a G [oojff]! ^ a is the identity. 
Hence the change of variables influences only the bilinear term in v. 
For brevity of notation, let 

V := Hl DQ {Si\T{a Q )) . 

On V we consider the norm || • ||y defined by \\v\\y := ||-Dt> H2, and the scalar product 
(v,w)v '■= (Dv,Dw), where || • || 2 and (•, •) denote the norm and, respectively, the 
scalar product in L 2 (ft) or L 2 (ft \ r((7 ));R 2 ), depending on the context. Let V 
denote its dual space and let (•, •) denote the duality pairing between V' and V. 
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For every t G [0, T], v G 7, and a G [(To, a] define 
F(t,v,a) := 

= <*j(<r)Dj(v + $(t))Di(v + ii>(t))dx + <T- g(t)(v + $(t))dH 1 . 

Jf2\r(o-o) ij - e | lj2 | Jd N Q 

Then the functional T can be also written as 

F(t,v,a):= (A(a)Dv\Dv)dx + 2 (Dip(t)\Dv) dx - 

Jn\T(ao) Jn\T(cro) 

-I g(t)vdH 1 + a + f \Di){t)\ 2 dx- f g(t)tp(t) dH 1 , 

or 

F(t,v,a) : = 

= (A(cr)Z^, + 2(£ty(t), £>u) - (</(t), i>Wn + o" + I WW 111 - (g{t),m)d N n 

where (•, -)q n ci denotes the scalar product in L 2 (<9Arf2, 7Y 1 ). Hence the elastic energy 
becomes jF eZ (t, v, a) := Fit, v,cr) — a, and there exist four positive constants Ajr, 
Ajp, and Mjf, independent of t and cr, such that for every t G [0,T] and every 
cr G [g-0,0 7 ] 

^(t,v,a) > Ajr|H|^-^ 
^ eZ (t,t;,(r) < A T \\v\\ 2 v + M T , 

for every v G 7. Indeed, this follows from the uniform ellipticity of the bilinear part 
and standard estimates (on Q + and 

2.7. Critical points of the energy. For every t G [0, T] the function F(t, •, •) : 7 x 
[(To, a] — > M. is twice Frechet partially differentiable with respect to (v, cr). The partial 
differential dyj-'it, v, cr) belongs to V, while the partial gradient grad„.F(i, v, a) is, by 
definition, the element of 7 given by 

(grad^(t, v, a),w) v = 2(A(a)Dv, Dw) + 2(ip(t),w) v - (g(t),w) dN n, 

for every w G 7. The partial differential d a T{t, v, a) is given by 

d a F(t, v, a) = (d a A(a)Dv, Dv) + 1 . 

For fixed i> G 7 and cr G [cr , a], we have that v, a) G W 1 ' 00 ^, T), with 

d t F(t, v, a) = 2(DiP(t), Dv + DiP(t)) - (g(t),v + iJ>(t)) dlfSl - (g(t),^(t)) dN n . 

Note that by the regularity assumptions on ip and g it follows also that the map 

(t, v, cr) (grad^(t, v, a), d c F{t, v, cr)) 

is continuous from (0, T) x 7 x (<7 , cr) into 7x1 

The second order partial differentials with respect to (v, cr) are given by 

((df v>a) F(t, v, cr)(w h n), (w 2 , r 2 ))) = (A(<t)£>iui, Dw 2 ) + (d a A(a)Dv, D Wi )t 2 + 

+ (d a A(a)Dv, Dw 2 )n + (d 2 aa A(a)Dv, Dv)r x r 2 , 

for every (wi, Tj) G 7 x R, i = 1, 2, where ((•, •)) denotes the duality product between 
7'xl and 7x1 

Since, for fixed t and a, the function v i— > jF(t, t>, cr) is strictly convex, it has a unique 
critical point i> tj0 ., and v t>a is a minimum point. Also the function u \— > £(t)(u, cr) is 
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strictly convex and its critical point is the unique minimum point u t , c G AD(ip(t), a) 
of u i — ► £(t)(u, a). The function u t>a satisfies 

2/ (Du t , a \Dw)dx= f g(t,x)wdH 1 Vw G H^Q \ F(a)) . 
Jn\r(a) Jd N n 

Proposition 2.1. For fixed t G [0, T] critical points ofj-'it,-,-) correspond to critical 
points of £(t) in the following sense: minimum points v t ^ G V of v 1— > JF(£, r,o) 
correspond by the change of variables to minimum points u t ^ G i4D(^(t),(r) 0/ -u 1— > 
£(t)(w,o-). Moreover, d a T{t,vt )(J ,a) = d a E{t,a), where E(t,a) := £(t)(u tt a, o~). 

Before giving the proof we discuss some properties of the minimizers u t>a . The 
following result provides a useful characterization of the "singular" part of the dis- 
placement u t>a near the tip 7(0") of the crack. For the proof we refer to [12], [13]. 

Proposition 2.2. Let a G [er , a] and u G H 1 (il\ T(a)) be such that 

Au G L 2 (Vt \ r(o-)) and d v u = on T(a). (2.6) 

T/ien there exists k G R satisfying 

~2 




u -K^r^sm-eH 2 (U\r(a)), (2.7) 

/or every [/ CC f2 open. In (2.7), r(x) := \x — 7(0") | and #(x) is the continuous 
function on U\T(a) which coincides with the counterclockwise oriented angle between 
7(0") and x — j(cr), and vanishes on the points of the form x — 7(0") + /17(a) /or /i > 
sufficiently small. 



The coefficient K\j2/n represents the stress intensity factor associated to the dis- 
placement u at the tip 7(0"). We shall use its following characterization. 

Proposition 2.3. Let a G [o" ,(x], u G iJ 1 (fi \ T(a)) satisfying (2.6), and let k be 
defined by (2.7). Then for every = (0i, 2 ) G C^°(Q; R 2 ) we have 



^(7(^)7(0-)= / \((D 1 u) 2 -(D 2 uf)(D l( f )l -D 2( f )2 ) + 
Jn 1 

+2D x u D 2 u{D x <\) 2 + D 2 <j>i) dx + 2 Aw(L>iw0i + D 2 u(p 2 ) dx. 

J in 



(2.8) 



Proof. For a complete proof see [1, Proposition 2.2]. The idea is to consider rj > 
such that 5(7(0-), 77) C f2, to integrate by parts: 

/ [((A^) 2 -p2«) 2 )Pl01-^20 2 )+ 2^^^(^102 + ^201)1^ ^ 

and to pass to the limit as rj — > using (2.6). □ 
Proposition 2.4. The function a 1— > E(t,a) is differentiable on [ao,o] and 

d a E(t,a) = l-Kl a , (2.10) 

where n t ^\J\ is the stress intensity factor associated to u t , a at 7(0). 

Proof. The proof follows the same arguments of the proof of [1, Theorem 3.3]. To 
compute the partial derivative d a E(t, a) we consider a diffeomorphism similar to $ CT 
and then use Proposition 2.3. □ 
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Similar computations have been recently done in [15] when the stored energy density 
W is a polyconvex function with = oo for every matrix A with det A < 0, and 

T is a segment. 

Remark 2.5. Fix t G ]0, T[. The map a i— > v to , a has the same regularity as 
a i— > A(a), hence, under the regularity assumptions we made on A(o~), it is of class 
C 2 ([<Jo, a]). Since in this case we are not interested in the dependence on t, let us 
simplify the notation and set v a := v t0ta - Then standard arguments for elliptic PDE's 
allow us to obtain that for every a* G [o"o,a] there exists v' aif G V as strong limit 
in V of the difference quotient V<J ~^J* , and the map a i— > v' c is continuous in the 
strong topology of V. The same arguments can be repeated to obtain that there 
exists v "» G V as strong limit in V of the difference quotient V "~_ V J* and that the map 
a I— > t>" is continuos with respect to the strong topology in V. Note that i/ and v" a 
solve the following equations 

(A(a)Dv' a , Dw) + (^((j)^, Dw) = VweV, 

(A(<t)Du", Dw) + 2(d a A(a)Dv' a , Dw) + (d 2 a A(a)Dv a , Dw) = Vu* G V , 

respectively. 

Proof of Proposition 2.1. It follows from the change of variables, Proposition 2.3, and 
Proposition 2.4. □ 

Remark 2.6. Fix to G ]0, T[. With the same notation as in Remark 2.5 i> CT := i>t , CT , 
note that the second order differential, d? v ^(to, v a , a), of T with respect to (v , <r) 
is strictly positive definite if and only if the second order derivative of the function 
o i— > jF(t , t> CT , <j) is strictly positive, when both exist. Moreover, by Proposition 2.1, 
this is equivalent to the fact that the second order derivative of a i— > E(t , a) is strictly 
positive. 

Indeed, as d v J r (t ,v (T ,o-) = 0, and a \— > v a is, by Remark 2.5, a C 2 -function, we 
have 

{d a d v F{t , v a , a),w) + (d^Ffo, v a , a)v' a , w) = MweV. 

Assume that 

< ^(d a J r (t ,v a ,a) + (d v F{to,v a ,o-),v' a )) = 

= d? a F(t , v a , a) + (dadyFito, v a , cr),v' a ) + (<9„<9 CT .F(£o, v a , o-),v' a ) + 
+(dl v F(t , v a , a)v' a , v' a ) + (d v F(t , v a , a),v") . 

Hence 

dl a F(t , v a , a) + (d v d a F(t , v a , a), v' a ) > , 

which implies that 

dl a F(t ,v a ,a) > (d^Ffavv^W^v'r) 

(recall that in our case {d^T^w) = (p v d a T , w)). 
Therefore 

((9( v ,a)H t o, v a ,a){w,T), (w,t))) = 

= d 2 a(J T{t^ v a , a)r 2 + 2(d a d v F{t , v a , o),w)t + (di^Ffa, v a , a)w, w) > 

> (d 2 vv F(tO: Va, CT)v' a , <>T 2 - 2(^(<o, V„ w)t + 

+ (d 2 v J r (t ,v a ,o-)w,w) > 0, 
which shows that d? v sJ 7 (to,v er , a) is strictly positive definite. 
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It is also easy to see that if d^^J 7 ^, v a , a) is strictly positive definite then the 
second order derivative of the function a h- > J 7 (to, v a , a) is strictly positive. 

3. Irreversible quasistatic evolution 

Given an initial crack length ctq > 0, and an initial value, uq, of the displacement, 
such that the initial configuration is in equilibrium, we want to study a quasistatic 
evolution of configurations (u, a) which starts from (uq, cr ). We are interested in the 
evolution until the crack length reaches the value o\. We cannot avoid the solution 
to have jumps (even at t — 0) to configurations with crack lengths larger than o~\, if 
this is the case, then the boundary data are not compatible with a progressive crack 
growth on the interval [<7o, <7i]. 

Definition 3.1. The irreversible quasistatic evolution problem consists in finding a 
left-continuous map t h- > (u(t), cr(t)), where a(t) represents the length of the crack up 
to time t, and the displacement u(t) belongs to AD(ijj(t), a(t)), which satisfies the 
following three conditions: 

(a) local unilateral stability: for every t 

E(t)(u(t),a(t)) <£(t)(u,a(t)) e AD(i{;(t),a(t)) (3.1) 

d a E(t,a(t))>0, (3.2) 

where E(t,a) is defined in Proposition 2.1; 

(b) irreversibility: the map t i— > er(£) is increasing; 

(c) energy inequality: for every < s < t we have 

£(t)(u(t),a(t)) < £(s)(u(s),a(s)) + 

+ [ (2 [ (Du(r)\Dip(T))dx - [ g^u^dH^dT, 

Js ^ io\r(a(T)) Ja N n ' 

In terms of the functional J 7 , the irreversible quasistatic evolution problem consists 
in finding a left-continuous function t h- > (f (t), cr(t)) which satisfies the following three 
conditions: 

(a^r) local unilateral stability: for every t 

grad^(t, a(t)) = 0, 

d„F(tMt),<r(t)) >0; 1 ' J 

(bjr) irreversibility: the map t i— > cr(t) is increasing; 
(cjr) energy inequality: for every < s < £ we have 

F(t, v(t), a(t)) < F(s, v(s), a(s)) + J d t F(r, v(r), a(r))dr . 

A solution, t i— > (f (t), cr(i)), to this problem is called an irreversible quasistatic evo- 
lution for T . 

Let us remark that, by the very construction of the functional J 7 , an evolution 
for T is well-defined only for cracks whose length is less than or equal to a. 

In terms of an irreversible quasistatic evolution t i— > (v(t),a(t)) associated to the 
functional J 7 , the Griffith's criterion can be expressed as: 

' a(t) > 

d a r(t,v(t),<r(t))>0 (3-4) 
d a F(t,v(t),<r{t))*(t) = 
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for a.e. t. Since the first two conditions are included in the definition of an irreversible 
quasistatic evolution, it remains to prove the last one. 

Proposition 3.2. Lett \— > (v(t),a(t)) be an irreversible quasistatic evolution for T . 
Then for a. e. t we have 

d a F(t,v(t),a(t))&(t) = 0. 

Proof. Since a is increasing, & exists at a.e. t. Fix t such that d-(t ) exists. As, given 
cr(t), the function v(t) is determined as the unique solution of gra.d v J-'(t,v,a(t)) = 0, 
the hypotheses we made on A(a) and on the data ip and g imply that w(to) exists, as 
strong limit in V of the difference quotient "(*)~^( to ) . 

As gradj,jF(t, v(t), o"(t)) = 0, from the energy inequality (cjr) we deduce that for 
a.e. t 

d a F(t,v(t),a(t))&(t) < 0. (3.5) 
Since a(t) > and d a F(t,v(t), a(t)) > 0, (3.5) implies the equality to be proved. 

□ 

Going back to the energy functional S, the Griffith's criterion now reads 

' a(t) > 

1 - K 2 (t) > (3.6) 
(1 - K 2 (t))a(t) = 



for a.e. t, where n(t)y | is the stress intensity factor associated to the displacement 

u(i) at the tip a(t) (see Proposition 2.2). Since by the change of variables we made, 
d a F(t,v(t),a(t)) = 1 - K 2 (t), the previous proposition shows that during an irre- 
versible quasistatic evolution the Griffith's criterion is satisfied. Note that this can 
be proved directly for S, following, for instance, the lines of [5, Theorem 6.1]. 

In the context of variational models for quasistatic crack propagation, the evolution 
of minimum energy configurations was studied (see, e.g. [4], [2], [8], [3]) and existence 
results were proved in a very general setting (see [3]). This kind of evolution is a 
solution to the following problem. 

Definition 3.3. The globally stable irreversible quasistatic evolution problem consists 
in finding a solution to the irreversible quasistatic evolution problem which satisfies 
the global stability condition: for every t 

£{t){u{t), a{t)) < £(t)(v, a) V<7 > a(t) Vv G AD(ijj(t), a) . 

During a globally stable irreversible quasistatic evolution the total energy is an 
absolutely continuous function of time and the energy inequality (c) becomes an 
equality. 

However, a solution to this problem is not completely satisfactory since, in order to 
get the global stability, we have to compare, at each time, the energy of a configuration 
with the energy of all admissible configurations with larger crack lengths. This is why 
we use another criterion of selection: among all irreversible quasistatic evolutions we 
choose the approximable ones, i.e. those that can be obtained as limits of solutions 
to a regularized evolution problem. 

In this paper we consider the regularized problem given by a modified e-gradient 
flow for the functional T . Since we are interested in an irreversible crack growth for 
a varying in the interval [aojfTx], we look for an increasing function a{t). Hence, we 
consider the positive part of the derivative of T with respect to a. Then, we modify 
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the evolution law for the crack length in such a way that it never reaches a. To this 
end we introduce a penalization factor A(cr) that can be any Lipschitz continuous 
function of a which is equal to one for a < a±, is strictly positive for ai < a < a, and 
is equal to zero for a = ~a. For instance, let 

;= fr- fry 

a — <j\ 

In such a way the evolution is the one given by the e-gradient flow, with the constraint 
that a is increasing, on the interval [ao, ai] that we are interested in, and it is modified 
by this artificial penalization term for a > a±, so that we do not consider it meaningful 
for a > G\. 

Definition 3.4. A function t h- > (v £ (t),o~ £ (t)) is called a solution to the initial value 
problem for the modified e-gradient flow for the functional T on [0, T] 

v £ = -gr&d v J 7 (t,v £ ,a £ ) 
ea £ = (-d a F(t, v £ , a £ )) + \(a £ ) , 
v £ (0) = u 
a £ (0) = (T , 

where A(cr) is given by (3.7), if t> £ G C 1 ([0,T]; V), a £ is a C 1 increasing function from 
[0, T] into [cr ,<7] and the first equation in (3.8) is satisfied in the following sense 

(ev £ ,w) v = -(grad v J 7 (t,v £ ,a £ ),w) v VweV V£e[0,T]. 

Note that (3.8) is a Cauchy problem for an ordinary differential equation in V x M.. 

Theorem 3.5. There exists a solution (v £ ,o~ £ ) to the initial value problem (3.8) with 
A(cr) given by (3.7), and the following energy estimate holds: for every s, t e [0,T] 
with s < t 

e f \\ve{r)\\ 2 v dr + e f dr + ^(f, < 

J s A((7 £ (rj) 

<T{s,v £ (s),a £ {s)) + J d t T{T,v £ (r),a £ (r))dT. (3.9) 

Proof. Taking into account the expressions of grad^jF and d u J r 1 the equations in (3.8) 
can be written as 

e{v £ ,w) v = -2{A{a £ )Dv £ ,Dw)-2{^{t),w)v + {g{t),w) dN n VweV 
ea £ = (-(d a A(a £ )Dv £ , Dv £ ) - 1)+A((7 e ) . 1 ' J 

Since the vector field defining the equation (3.10) depends on t only through the 
boundary data ip and g, it is Lipschitz continuous in t. Moreover, for fixed t, standard 
estimates show that it is Lipschitz continuous and bounded on the bounded subsets 
of V x R. Hence classical results on ODE's (see, e.g. [7]) give the local existence and 
the uniqueness of the solution. Since there exist a G C([0, T]) and (3 > such that 

(-grad^f , v, <t),v) v + a{-8 a T{t, v, a)) + X(a) < a(t)(\\v\\ 2 v + a 2 ) + (3 



for every (v, a) G V x R, the solution is defined on the whole interval [0, T]. 
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The function t h- > ^(t, v e (t), cr £ (t)) is then Lipschitz continuous on [0,T] with de- 
rivative given for a.e. t G [0, T] by 

j t F(t,v £ (t),a £ (t)) = d t F(t,v £ (t),a £ (t)) + (gmd v F(t,v £ (t),a £ (t)),v £ (t)) v + 

+d IJ F(t,v £ (t),a £ (t))a £ (t). 

Taking into account the equations satisfied by v £ and a £ , for every s, t G [0, T] with 
s < t we have 

u e (i), <r e (i)) - .F(s, u e (s), cr £ (s)) = 

= jT (&F(r, Ve (r),<r e (r)) -e||t> e (r)||^ ~ ^J^^yj > 

which implies (3.9). □ 

Remark 3.6. Let t i— > (v £ (t) , o~ £ (t)) be a solution to (3.8). Assume ||u e (i)||y < M for 
some positive constant M independent of t and e. By (2.3), 

ea £ (t) < (MM 2 + l)X{a £ {t)) < C(a - a £ (t)) + , 

for some constant C > 0. By classical results on differential inequalities (see, e.g. [14, 
Theorem 1.6.1]) it follows that for every t G [0, T] 

a £ (t) <a- e - ct / £ (a-a ), 

hence a £ never reaches a. 

Note that, since the evolution is constrained to cracks with lengths less than or 
equal to a, Griffith's criterion is meaningful in this setting only until the length a is 
reached. As the penalization factor X(a) is strictly positive for a < a, we may replace 
(3.4) by 

' a{t) > 

d a F(t,v(t),a(t))X(a(t)) > 
d <T F(t,v(t),a(t))&(t) = 0. 

for a.e. t G [0, T\. Therefore, also the second line in the local stability condition (ajr) 
may be replaced by d a J-(t, v(t), a(t))X(a(t)) > 0. 
We introduce now the following notion of evolution. 

Definition 3.7. The approximable irreversible quasistatic evolution problem on the 
interval [0, T] with initial data (u ,a ) consists in finding a left-continuous map t i— > 
(v(t),a(t)) from [0,T] into V x M which satisfies the following conditions: 

(a^-) for every t G [0,T] 

grad^(M(t),a(t)) = 
d a F(t,v(t),<j(t))X(o-(t)) > 0; 

(bjr) the map t \— > <r(£) is increasing; 
(cjr) for every < s < t < T 

.F(t,u(t),<7(f)) <^(s,u(s),<t(s))+ / ^(r,u(r),<r(r))dr; 
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(d^-) the pair (v(t), cr(t)) is the limit of a suitable solution (v e (t), cr e (t)) of the mod- 
ified e-gradient flow for T with initial conditions v £ (0) = u and cr £ (0) = a , 
in the sense that for a.e. t 

a £ (t)^a(t), , 3U , 
v £ (t) — > f (£) strongly in V. ^ ' ' 

A solution t t— > cr(t)) to this problem is called an approximable quasistatic evo- 
lution for JF. 

We are now in a position to state the main result of this paper. 

Theorem 3.8. There exists a solution t i— > (v(t),cr(t)) to the approximable irre- 
versible quasistatic evolution problem with initial condition (uo,a ) on [0, T]. 

Remark 3.9. The fact that an approximable quasistatic evolution starts from (uq, o~o) 
means only that for every e > 0, v £ (0) = uq and cr £ (0) = a . We may always 
set (f(0),a(0)) := (1*0,00), but in general v and a are not continuous in t — 0. 
The only case in which (u , a ) is the initial value for the evolution in a "classical" 
sense, is when (u ,a ) is the absolute minimum point of JF(0,-,-). Indeed, in this 
case, by semicontinuity and by the energy inequality (cjr), it is easy to see that 
t 1— > J-(t, v(t), o~(t)) is continuous in t — 0. 

Proof of Theorem 3.8. For e > let (v £ , a £ ) be the solution of the modified e-gradient 
flow with initial data (u ,a ). Let t G [0,T]. The estimates we have on together 
with (3.9) between s = and t imply 

AHMOIIv<^ + ^(0,«o,<to)+ / (a(r)||u e (r)||^ + 6(r))dr 

./o 

for some functions a, 5 G L°°(0,T) which depend only on the data ip and Then, 
by Gronwall's Lemma, there exists a positive constant C > independent of t and £, 
whose value may change from line to line, such that 

\\v £ (t)\\ v <C VtG[0,T]. (3.12) 

By (3.9) we now get 

£ \\ve\\ 2 L 2 i0 , T]V) < C (3.13) 

e||*.|li» ( o,r) < C ( 3 - 14 ) 
Let e — > 0. By Helly's Theorem, there exists a subsequence, still denoted by £, and 
an increasing function a: [0, T] — > [(To,?] such that 

<7 e (i) -> <r(i) for every t G [0, T] . 

The estimate (3.12) implies that there exists a function v G L 2 (0,T; V) such that 

v £ ^v weakly in L 2 (0,T; V) , 

while, by (3.13), 

ev £ -> strongly in L 2 (0, T; V) . 

Hence 

e(u e (t),u;)y = (-grad w J r (*,u e (t),<r e (t)),«;)v -> (-grad^f, v(t), a(t)), w) v = 0, 
for every w G V and for a.e. t G [0,T]. It follows that 

rT pT 

(A{a £ {t))Dv £ {t),Dv £ (t))dt^ / (A(a{t))Dv(t),Dv(t))dt, 



I 
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which gives the strong convergence in V of v e (t) to v(t) for a.e. t G [0, T]. 

By (3.14), ea £ (t) — > for a.e. t G [0, T\. Taking into account the equation satisfied 
by a £ , we obtain that {-d a J r {t,v(t),a(t))) + X(a(t)) = for a.e. t G [0,T]. 

When passing to the limit in (3.9), we neglect the terms containing the norms of 
the time derivatives of v e and a e , and thus get that for a.e. s, t G [0,T] with s < £ 

^(t, (j(0) < ^(s, v(s), a(s)) + jf <9^(r, v(r), a(r))dr . (3.15) 

(By semicontinuity the estimate holds true for every t G [0, T].) 

Since a is increasing, for every t G [0, T] there exists the limit o~ G (t) := lim s ^_ cr(s). 
Let v e (t) be the unique solution to grad„.F(t, t>, <j e ) = 0. Then v(s) — > t> e (t) strongly 
in V as s — > t— , er(£) = <r e (t) and f(t) = t> e (t) for a.e. £ G [0,T]. By construction, 
the map t \— > (f e (t), <r e (t)) is left-continuous from [0,T] into V x [<7 ,(7]. Moreover, 
9 (T J r (t,w e (t),a e (t))A(cr e (t)) > for every * G [0,T]. Let s, i G [0,T] with s < f, and 
let s n — > s— , £ n — > £— be such that (3.15) holds for s n and t n . Passing to the limit in 
(3.15) asn^ +oo we obtain 

F(t, v Q (t), a Q (t)) < F(s, v Q (s), a e (s)) + J* d t T(r, v Q (r), a Q (r))dr , 

so that we conclude that (v e , <r e ) is an approximable quasistatic evolution for T on 
[0, T] which starts from (u ,a ). □ 

Remark 3.10. From (3.15) we deduce that if t G [0, T\ is a discontinuity point of 
F(t,v(t),a(t)) then 

lim JF(£, o-(t)) < JF(£, tr(t)) . 

Indeed, note that at every time £ the function t h- > cr(£) has a right limit. Let 
<7®(t) := lim t _j + cr(t), and let w e (t) be the solution to grad^J 7 ^, <7®(i~)) = 0. By 
the regularity assumptions made on the data, we have that v(t) converges to v®(t) 
strongly in V, and hence, using (3.15), we obtain 

lim F(t,v(t),a(t)) =F(t,v®(t),a®(t)) < F(t,v(t),a(t)) . 



4. Quasistatic evolution and the Implicit Function Theorem 

In this section we show that, under suitable regularity assumptions, the solution to 
the modified e-gradient flow converges to the continuous solution for the quasistatic 
evolution problem given by the Implicit Function Theorem. 

Theorem 4.1. Assume that in (t°, a ) G [0,T[ x [cr , <Ji[ the following conditions are 
satisfied 

d a E(t°,a°)=0 (4.1) 
d 2 a E(t°,a°)>0. (4.2) 

Then there exists a time interval [t°,t l ] and a unique Lipschitz continuous function 
a : [tV 1 ] -> [<t°,<Ti] such that 



d a E(t,a o (t)) = Vtelt^t 1 ] 
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Moreover, if (v £ , a £ ) is the solution to the modified e-gradient flow and the following 
two conditions are satisfied: 

& £ (t)>0 VtGftV 1 ] 

then a £ {t) -> a°(t) and E(t,a £ (t)) -> E(t,a°(t)) for every t G [t ,* 1 ]. 

The first part of the theorem follows from the Implicit Function Theorem. As for the 
second part, let us remark that even if there are not at the moment general theorems 
guaranteeing the strict monotonicity of a £ during the approximation process, in many 
cases this will follow, for a suitable choice of the boundary data, from a symmetry 
argument. 

We now prove the second part of the theorem in an equivalent form for the func- 
tional T . Indeed, since d a E(t,a) = d a J-'(t,Vt i cr,o-) (see Proposition 2.1), if the sec- 
ond order derivative d^E(t°,a ) > 0, then also ^d a !F{t , v t o jtT o, a ) > 0, and this 
last condition is equivalent to the fact that the second order partial differential 
d? V(J \F{t°, v t o jtT o, a ) is strictly positive definite (see Remark 2.6). 

Theorem 4.2. Assume that in (t°,v°,a ) G [0,T[ x V x [<7o,<7i[ the following condi- 
tions are satisfied 

grad^(tW°) = 0, 
d a F(t o ,v o ,a o ) = 0, 

and the second order differential, d? a \J~(t°, v°, a ), of J- with respect to (v,o~) is 
strictly positive definite, i.e. there exists a > such that 

((^^^OKr),^,^))^^!^!!^^! 2 ) VweV VrGl. (4.3) 

Let (v £ ,a £ ) be the solution of the modified e-gradient flow for T given by Theo- 
rem 3.5 and assume that 

v £ (t°) — > v° strongly in V and 
a £ {t°)^a° ase^O. 

Then there exist a time interval [t ,^ 1 ] and a unique Lipschitz continuous function 
(v°, a ) : [t°, i^Vx [a , a ± ] such that 

(gr8,d v F(t,v (t),o-°(t))=0 

\d r(ty(t),o°(t)) = o 

for every t G tf 1 ]. Assume that & e (t) > and a £ (t) < U\ for every t G [i ,^ 1 ]. Then 
v £ (t) — > v°(t) strongly in V and cr £ (t) — > a°(t) for every t G [t ,^ 1 ]. 

Proof. By our assumptions on the data, d} va -.J-{t,v,o) (see Subsection 2.7) is con- 
tinuous with respect to (t,v,o~) G [0,T] x V x [o"o,(f]. Moreover, the function 
t i— > (9 t grad^^ r (t, v, a) belongs to L°°(0,T;V), while dtd^Fit, v, a) = 0. By the Im- 
plicit Function Theorem (see, e.g., [16]) applied in (t°,v°, a ) to 

grad„jF(t, v, a) — 
d a T{t,v,o) = 0, 
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it follows that there exist a time interval [t ,^ 1 ] and a unique Lipschitz continuous 
function (v°,a°): [t ^ 1 ] -> Vx[<t°,<Ti) such that 

'grad^(M°(t),a°(t)) = 

^(M°(t),a (t)) = 1 ' J 

for every t G [i , i 1 ]. By a compactness argument, changing eventually the value of a, 
we may assume that there exist a > and r > such that for every t G [t ,^ 1 ], for 
every v G B r (v°(t)) C V, and for every a G (<x°(£) — r, cr°(i) + r) 

((,9 ( 2 , )(T) ^(t,t;,a)( W ,r),( W ,r)))>a(||«;|| 2 +|r| 2 ) V^V VrGl. (4.5) 

Restricting eventually the time interval, we have <J°{t) + r < o\ for every 

Let < r' < r be a number that we shall choose later. For every e > small 
enough we have ||v £ (t°) — v°\\v < r ' and \cr e (t°) — o°\ < r' . By continuity, there 
exists a time interval, depending on e, on which these inequalities hold. Let r £ be the 
largest time such that for t < r £ , \\v £ (t) — v°(t)\\y < r' and \cr e (t) — cr°(t)\ < r' . Then 
<J £ (t) < (Ti, hence X(a e (t)) — 1 for t < r £ . 

We want to prove that t £ — t 1 . Assume by contradiction that r £ < t 1 . Taking 
v £ (t) —v°(t) as test function in the equation satisfied by v £ , multiplying by a £ (t)—a°(t) 
the equation satisfied by a £ , and taking also into account (4.4), we obtain 

= -(grad^(t, v £ (t), <r e (t)) - grad^(t, v°(t), a°(t)),v £ (t) - v°(t)) v + 
+ ( - d a F(t,v £ (t),a £ (t)) + d^(t,v°(t),a°(t)))(a £ (t) - a°(t)) - 
-e(v (t) : v £ (t)-v (t)) v -ea (t)(a £ (t)-a°(t)). 

Setting 

Ut) := \\v e {t)-v°mv + Mt)-o°(t)\ 2 , 
from (4.5) it follows that 

|4W < -a( £ (t)-e(v (t),v £ (t)-v (t)) v -ea (t)(a £ (t)-a°(t)) < 

< -<*Ut) + \ ll*°(*)llv + §IM0 - «°(*)llv + f l*°(*)l 2 + §M*) - ^°(0I 2 < 

<(- a + £)C £ (0 + |/3 WG[t°,r £ ), (4.6) 

where (5 is an upper bound for + |<7°(t)| 2 on [t ,^ 1 ]. 

Hence 

Ut) < (Ut°) - ^K- +1)(t -< 0) + ^7 Vt G [f,r e ) . (4.7) 

Therefore, choosing now r' small enough, from (4.7) we get that also ||t> £ (r £ ) — 
^ ( r e)||v < r an d \<j £ (t £ ) — o"°(r £ )| < r. By continuity, these inequalities hold also 
for some t > r £ , which contradicts the maximality of r £ , and so we deduce that 
r £ = t 1 . Hence (4.7) holds for every t G [t ,^ 1 ]. Passing to the limit in (4.7) we get 
the conclusion. □ 

By the change of variables that defines the functional J 7 , and by the uniqueness 
of the regular evolution given by the Implicit Function Theorem, it follows that the 
regular evolution in Theorem 4.2 corresponds to the one in Theorem 4.1. 
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Proof of Theorem continued. Let (v £ , a £ ) be the solution to the modified e-gradient 
flow. By Theorem 4.2, a £ (t) — > <r°(t) and v £ (t) — > v°(t) strongly in V for every 
t G i 1 ]. Since the function v i— > grad„jF(t, t>, a) is continuous from V to V with 
respect to the strong topology, it follows that 

grad^f , v £ (t), a £ (t)) -> grad^f , <r°(t)) = . 

Let v e {t) be the element of V associated to u t ^i t \ by the change of variables. As 
grad 1 ,jF(t, v e (t), <r £ (t)) = we deduce that v £ (t) — v £ (t) — > strongly in V. This 
implies that 

.F(t,U e (t), a £ (t)) - F(t, v e (t), a £ (t)) -> , 

On the other hand, 

.F(*,u e (t),<7 e (f)) -.F(*y(*),<7 (*)) = E(t,a°(t)) 
F(t,v £ (t),a £ (t)) = E(t,<j £ (t)) , 

so that we conclude that £7(t, <r e (t)) -> 0-°(*)) for every t G [t°, i 1 ]. □ 
5. MONOTONICALLY INCREASING LOADINGS 

In this section we consider the setting proposed by Francfort and Marigo in [9] and 
compare the evolution defined therein with a solution to the irreversible quasistatic 
evolution problem. 

Assume ip(t) := tipo, with tp G and g(t) = 0, and define 

E{a) := min{||DM||2 : u G AD(ijj ,a)}. 

Since H\Q \ T(a')) C H\Q \ T(a")) for a' < a", we have that E(a') > E{a"), so 
that the function a \— > is decreasing. 

Following Definition 4.13 in [9], we define a crack trajectory t i— > a FM (t) by the 
following three properties: 

(i) £ i— > a FM (t) is increasing; 

(ii) t 2 E(a FM (t)) + o> M (i) < t 2 £(cr) + a for every tr > 0> M (*); 

(iii) t 2 E(a FM (t)) + a FM (t) < t 2 £(<7 FM (s)) + a FM (s) for every s < t. 

The following result shows that if a — > -E'(cr) is concave in some subinterval of (do,?) 
then t i— > a FM (t) is discontinuous. 

Proposition 5.1. Let £ i— > a FM (t) be a crack trajectory which satisfies properties (%)- 
(iii) above. If there exists a subinterval (a,b) C (o~q,(t) where o~ \— > £'(cr) zs concave, 
then o FM (t) has some discontinuity points. 

Proof. Let £o > be such that o~ FM (to) < a. We first prove that there exists t > t 
such that a FM (t) > a. Indeed, assume by contradiction that a FM {t) < a for every 
t > to- Then conditions (i), (ii) and the fact that a h- > -E'(cr) is decreasing imply the 
following inequalities: 

t 2 E(b) + b > t 2 E(a FM (t)) + a FM (£) > t 2 E(a) + <r FM (£) > t 2 E(a) + a FM (t ). 

In particular, we deduce that t 2 < b(E(a) — which, up to considering T large 

enough, represents a contradiction. 

If cr FM (t) 7^ a for every £ G [0, T] then o> M is discontinuous and the proof is 
concluded. Otherwise, let £ be the first time such that cr FM (t) = a. We claim that 

0>m(*) = a for every £<£<£*, (5.1) 
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where 

^ " V 1 E(a)~-E(by (5 ' 2) 

Indeed, fix t e (M*) and assume by contradiction that a FM (t) £ ]a,b]. Then there 
exists a e]0, 1] such that a FM (t) = aa+ (1 — By condition (ii) and the concavity 
of cr ^ E(a) we have 

t 2 E(b) + b> t 2 E(a FM (t)) + a FM (t) > t 2 aE(a) + t 2 (l - a)E(b) + aa + (1 - a)fe, 

that is 

t 2 £(a) + a < t 2 £(6) + 6. (5.3) 

Therefore 

t 2 E(a FM (t)) + a FM (t) > a(t 2 E(a) + a) + (1 - a)(t 2 £(6) + 6) > t 2 £(a) + a. (5.4) 

Since (5.4) is in contradiction with condition (iii), we deduce that cr FM (t) = a. 
Consider now the case t = t*. Formula (5.3) becomes the identity 

(t*) 2 E(a) + a = (t*) 2 E(b) + b. 

Assume there exists a e]0, 1] such that a FM (t*) = aa + (1 — a)b; then, arguing as 
before, we obtain that 

(t*) 2 E(a FM (t*)) + <j FM (t*) > (t*) 2 E(a) + a, 

which, by conditions (ii) and (iii), implies that a FM (t*) = a. 
To conclude, we prove that 

a FM (t) >b for t > t*. (5.5) 

Indeed, let us fix t > t* and assume by contradiction that <7 FM (t) < b. Then there 
exists a g]0, 1] such that a FM {t) = aa + (l — a)b, and this fact together with condition 
(ii) implies that 

2 b - a FM (t) a(b-a) = 2 

" E(a FM (t)) - E(b) ~ aE(a) + (1 - a)E(b) - E(b) 1 } ' K ' } 

a contradiction. This fact concludes the proof, since we have shown that for t < t* 
Cfm(0 = a, while a FM (t) > b for t > t*. □ 

Let («(•), cr(-)) be an irreversible quasistatic evolution. Recalling that u(t) is the 
minimum point of ||-Dm|| 2 on AD(tip ,a(t)), we have that = t 2 E(a(t)). 

We may now express conditions (a), (b) and (c) of Definition 3.1 of an irreversible 
quasistatic evolution, in terms of cr(t), and, in the case of this particular choice of the 
data, we obtain: 

(a') 1 + t 2 E'(a(t)) > for every t > 0; 
(6') the map t i— > a(t) is increasing; 

(c') t 2 E(a(t)) + a(t) < s 2 E(a(s)) + a(s) + 2 tE{o{t)) dr, for every < s < t, 

J s 

where E'(a(t)) denotes the derivative of E with respect to a computed at a it). 
Since E{a{r)) < E(a(s)) for r G [s,t], condition (c') implies condition (iii). 

Remark 5.2. Let t h- > a(t) be a left-continuous map on [0, T] which satisfies condition 
(c') and define 

a G {t) := km sup— — . 

s^t~ t — S 
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Then 

(l + t 2 E'(a(t)))a e (t) <0 (5.7) 
for every t G [0, Tj. Indeed, let t k / t be such that 

hm a{t) ~ a M = & e {t y 

k^oo t — t k 

Then condition (c') between t k and t can be written as 

{t 2 - t 2 k )E(a(t)) + t 2 k (E(a(t)) - E(a(t k ))) + a(t) - a(t k ) < 2 f rE(a(r)) dr, 
and (5.7) follows dividing by t — t k and letting k — > +oo. 

Remark 5.3. Let t i— > a(t) be a left-continuous map on [0, T] which satisfies condi- 
tions (a'), (6'), and (c'), and let t > be such that <? e (t) > 0. Then, by Remark 5.2 
and conditions (a') and (6'), it follows that 

E'(a(t)) = ^-E(a)l =a(t) = -^, 

which implies that a(t) does not belong to the concavity intervals of E(a), since 
t i— > cr(t) is increasing, and t i— > E'(a{t)) would be decreasing, while the right-hand 
side is increasing. More precisely, if there exists an interval (a,b) C [<7 ,cx] such that 
a i— > E'(cr) is strictly decreasing on (a, 6) and there exists t > such that cx(to) > 
(or o" e (t ) > 0) and o"(t ) G (a, 6) then we reach a contradiction. Indeed, let t > t 
be such that a{t) G (a,b). By (6'), <x(t) > cr(t ), and by (a') and our assumption on 
E'(a), we get 

1 _ 1 1 



--<E'(a(t))<E'(a(t )) = -^<- t2 . 



a contradiction. 



In order to specify better the monotonicity needed in the above remarks we intro- 
duce the following notion. We say that to is a local left- constancy point for a if there 
exists e > such that a is constant on the interval [t — e,t ]. 

Proposition 5.4. Let a: [0, T] — > [o"0)(x[ &e a left- continuous map which satisfies 
conditions (a'), (b'), and (c') 7 and let t > 0. // 

(1) to is not a local left- constancy point for a and 

(2) there exists (a, b) C [<7o,ct[ s,uc ^ -^'(cO ^ s strictly decreasing on (a,b) 
then a(t ) ^ (a, 6). 

Proof. If to is n °t a local left-constancy point for a, then, given e > 0, there are 
t\,t 2 G [to — e,to] such that cr(t^) ^ cr{t 2 £ ). Therefore, there exists t £ G [to — £,to[ such 
that d e (t £ ) > 0. Then (5.7) together with (a') imply that 1 + t 2 E'(a(t £ )) = 0. By 
Remark 5.3, cr(t e ) (a, b) and we conclude by passing to the limit as e — > (since a 
is left-continuous). □ 

Proposition 5.5. Let a : [0,T] — > [cr ,a[ be a left- continuous map which satisfies 
conditions (a'), {b'), and (c'). Assume that E(a) is convex on (a, b) C [oojO 7 ]- T/ien 
er(t) continuous at every t with o~(t) G (a,b). 

Proof. Assume by contradiction that a(t) < o"(t + ). Then condition (c') and condition 
(a') imply 

E(a(t+) - E{a(t))) 1 

a(t+) - a(t) " "i* " E 
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a contradiction. 



□ 



6. Concavity and convexity intervals for the energy functional 

We consider the energy functional 

a ^ E(a) : = min{||Dzi||2 : u G AD(ip,a)}, 

and construct an explicit example of Q and ip for which E(a) is concave on some 
subinterval. Let -B_2 denote the ball of radius 1 centred in (—2,0), let B2 denote the 
ball of radius 1 centred in (2,0), and let F : = [—3,3] x {0}. 



B-2 




B 2 

( /v\ 




\ -2 J 




\ 2 / 






Figure 1. The set fL 



B_ 2 UT £ UB 2 



For e > let 

T £ := ] — 2 + cose, 2 — cose[ x ] — sine, sin e[, Vl £ 

Further, for every a G [—3, 3] let 

T(a) := [-3, a] x {0}. 

Let (p, 9) and (p, 9) be polar coordinates around (—2, 0) and (2, 0), respectively, where 
the functions 9 and 9 are chosen, as in Proposition 2.2, such that 9(xi ) x 2 ) — > — vr if 
x 2 — > 0— and x\ < —2, 9(xi,x 2 ) — > 7r if rr 2 — ► 0+ and xi < —2, and, analogously, 
X2) — > — vr if rr 2 — * 0— and X\ < 2, x 2 ) — * 7r if X2 — > 0+ and xi < 2. 
On <9f2 e we define the boundary data ip £ as follows: 

on (dB_ 2 r)dn £ ) \ T(a), 

on (dB 2 ndn £ )\T(a), 



sin 
sin 



6_M 

2 

2 



^ e (x) := < 



sin § 



— sin 
sin 

— sin 



£ 

2 

£ + 



2 1 2-cose 
£ _l_ ^1 
2 



cos I — sin |) 



: ( sin I — cos 6 



on ] — 2 + cose, 0[ x {sine}, 
on ] — 2 + cose, 0[ x {— sine}, 
on [0, 2 — cos e[ x {sin e}, 
on [0, 2 — cos e[ x {— sin e}. 



(6.1) 



2-cose V 2 2) 

For every a g] — 3, 3[ , let w £ (cr) G B x {£l £ \ T(cr)) be the solution of the problem: 
E (a) := mm / \Du\-(lv : 11 G LD(V .a) ^. (6.2) 



{ / |-Dm| 2 dx : u £ AD{i[) e , u)\. 

1 JQ £ \r(a) > 



Our aim is to prove that for e sufficiently small there exists a subinterval [a, b] of 
[—2,2] such that E £ (a) is concave on [a, b]. 

As (7 1— > E £ (cr) is a C 2 -function, in order to prove that E £ (<j) cannot be convex on 
the whole interval [—2, 2], it is enough to show that the following three conditions are 
satisfied: 

(a) limsup e ^ 0+ E £ (2) is finite; 
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sin ■ 



sin ■: 



2 — COS € 



(cos | - sin | ) 



sin : 



sin ; 



\ r(<7) 



sin | 



- sin l+2^o77( sill i- C0S |) 



Figure 2. The boundary datum ip £ . 



(b) liminf e ^ + E e {-2) = oo; 

(c) limsup e ^ + E' e (—2) is finite; 

where we denote by ' the first derivative with respect to a. 

In order to prove condition (a) we construct an admissible function u £ for E £ (2) 
whose energy, ||Dii e |||, is bounded uniformly with respect to e. We define the open 
sets BZ 2 and BZ 2 by 

B -2 = {(xi,x 2 ) e B_2 ■ x 2 > 0} 
B ~2 = {(^1,^2) e B_ 2 : x 2 < 0}. 

Let v + be the solution to the following problem: 

Au = on Bt 2 ; 

u(x) = sin ^ on dB+ 2 n 95_ 2 ; 
9„u = on ] - 3, -1[ x {0}. 

Then the function v~{x\,x 2 ) '■— —v + (xi,—x 2 ) solves the analogue problem on BZ 2 - 
Let u £ be the function which coincides with the harmonic functions that satisfy the 
boundary conditions on BZ 2 , on BZ 2 , and on B 2 , respectively, that is, u £ := v + on 

B~Z 2 i u £ := v~ on BZ 2 , and u £ := sin | on B 2 . On T £ \(B 2 LiB_ 2 ) we define m £ in the 
following way: on the horizontal line x 2 = sin#, with 9 G [—£,£], we set u £ (xi,x 2 ) : = 
sin I for x\ e] — 2 + cos#, 0] and then interpolate linearly with the boundary data on 
dB 2 nT £ : u £ {x u x 2 ) := sin | + ^^(cos f - sin |) for Xi G [0, 2 - cos 6[ , if < 9 < e, 
and u £ (xi,x 2 ) := — sin | + 2 _^ sg (— sin | — cos f ) for xi G [0, 2 — cos^f , if — e < 9 < 0. 
It is easy to check that u £ G AD(ip £ ,2) and that Du £ is bounded in L 2 (f2 £ \ r;M 2 ) 
uniformly with respect to e. This implies that 

limsupi? £ (2) < limsup / \Du £ \ 2 dx < +00, 

and condition (a) is satisfied. 

We continue by proving condition (b), i.e., E £ {—2) tends to infinity as e goes to 
zero. Let us first consider the model problem 

min { / \Du\ 2 dx : u > - ondiR £ ,u<—- on d 2 R £ } (6.3) 
JRe 2 2 

where 



R £ :=]0,l[x] -e,e[, d 1 R £ := [0,1] x {e}, d 2 R £ := [0, 1] x { - e). (6.4) 
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It is easy to see that problem (6.3) admits a solution and that it is equivalent to 

[11 

min { / \Du\ 2 dx : u = — ondiR £ ,u = —- on d 2 R £ }, (6.5) 

J Re 2 2 

which admits the affine solution u a (xi,X2) '■— 2 \ x 2 f° r ever y x = ( x i, x 2) £ -Re- 
Going back to the domain Q £ , let us consider the same problem with different 
constants: the rectangle R £ is defined now by 

R £ := }A £ , 2 — cose[ x ] — sine, sine[ C T £ , 

where A £ is a positive constant such that ip £ (x) > \ on d\R £ := [A £ , 2 — cose] x {sine}, 
(and i> £ (x) < —\ on d 2 R £ :— [A £ ,2 — cose] x {—sine}), when e is sufficiently small. 
Then 

E e (-2) = [ \Du £ {-2)\ 2 dx> [ \Du £ {-2)\ 2 dx> ( \Du a \ 2 dx . 

Jn e \r(-2) JR e JR e 

Since J R \Du a \ 2 dx — > +oo as e — > 0, condition (b) is proved. 

It remains to show that condition (c) is satisfied, i.e., that the first derivative of 
a i— > E £ (o) at a = — 2 is bounded as e goes to zero. Since 

E>(a) = -K 2 £ (a), (6.6) 



see, e.g. [13, Theorem 6.4.1], where K £ (a)^J | is the stress intensity factor associated 

to u e (a) at the tip (ex, 0), see Proposition 2.2, it is enough to show that n £ (cr) remains 
bounded when, for instance, — | < a < — |. 

For a E [—5/2, —3/2], let v(o) be the solution of the following problem: 

f 9 
min{ / \Du\ 2 dx : u eH l (B^ 2 \V {a)), u = sin- on 9fi_ 2 \ T(a)}. (6.7) 

Let us extend v(a) to R x [—1, 1] constantly on the horizontal lines and denote now 
by v{a) this extension. 
We claim that 

u £ (a) -> v{a) strongly in H\B_ 2 \ r(<r)). (6.8) 

Assuming the claim true, we now use the following characterization of k £ (see 
Proposition 2.3): 



n 2 e {a) 



[ [((Aiif - (D 2 u £ ) 2 )D lV + 2D lU £ D 2 u £ D 2 J dx (6.9) 



with (/? G C c 1 (-B_ 2 ) such that </?(cr, 0) = I. By (6.8) and the definition of v(a), we 
can pass to the limit in the right-hand side as e — > + and define in such a way the 
quantity: 

n 2 (a) := [ IdD^ia)) 2 - {D 2 v{a)) 2 )D lV + 2D 1 v(cr)D 2 v(cr)D 2 J dx. (6.10) 

JB- 2 \r(a) 1 J 

Therefore, by (6.6), 

5 3 

limsup E' £ (a) = —k 2 (<j) for every < a < — . (6-11) 

£ ^o+ 2 2 

As, by (6.10), k{o) is bounded, formula (6.11) concludes the proof of condition (c). 



24 RODICA TOADER AND CHIARA ZANINI 

Proof of the claim. Let fi £ := T £ U B 2 and let w £ be the solution of the following 
problem: 

min{J \Du\ 2 dx:ueH\n £ ),u = ij £ on <9fi £ H <9fi e } . (6.12) 

We consider a cut-off function tp e C°°(R) such that < tp < 1, = 1 for 

X! < — |, and = for X\ > — |. Then the function ( := tpv(a) + (1 — tp)w £ 

belongs to AD(ip e ,a) and 

E £ (a)= [ \Du £ (a)\ 2 dx< [ \D(\ 2 dx. (6.13) 
By convexity, we have 

/ \D(\ 2 dx< [ \Dv(a)\ 2 dx+ I \Dw £ \ 2 dx + [ \Dv{a)\ 2 dx+ 

Jn e \r(a) JB_ 2 \r(o-) Jn e Jt e 

+ / (2Dtp(tpDv(a) + (1 - tp)Dw £ )(v(a) - w e ) + \Dp\ 2 { v(a) — w £ ) 2 ) dx. 

JT £ n(suppDip) 

(6.14) 



Now 



and, for any r] > 0, 



lim / \Dv(a)\ 2 dx = 0, 
lim / \Dtp\ 2 (v(a) - w £ ) 2 dx = 0, 

: ^ 0+ J T E C\{suppDu>) 



(6.15) 



/ 2Dip(ipDv(a) + (1 - y?)£>w £ ) (u(<r) - w e ) dx < 

J T £ n(suppDifi) 

<2 D(p(pDv(a)(v(a) - w £ ) dx+ (6.16) 

JT £ n(suppDip) 

+- \Dp\ 2 \v(a) - w £ \ 2 dx + 7] I \Dw £ \ 2 (l - p) 2 dx . 

V JT £ n(suppD(p) JT £ n(suppD(p) 

Since the first two terms in the right-hand side tend to zero, it remains to prove that 

lim / \Dw £ \ 2 dx = 0. (6.17) 

e^0+ J T £ n(suppD<p) 

As in the proof of condition (b), we consider first a model problem. Similarly to (6.4), 
we now set 

R £ :=]-l,0[x ]-e,e[, d ± R £ := [-1,0] x {e}, d 2 R £ := [-1, 0] x {-£>, 
and define /i £ as the solution to the following problem: 

'Ah e = on R £ , 
h £ = \ on (9ii? £ , 

h e = -f on <9 2 -R £ , 



(6.18) 



IMloo < 1- 



We claim that 



lim / \Dh £ \ 2 dx = 0, (6.19) 



Re 



where 

i 4 1 r , 
R £ :=] - -,--[x]-e,e[c R £ 
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Indeed, note that the function z e (xi,x 2 ) '■= \x 2 solves (6.18) (for e < 1). By a 
Cacciopoli type estimate we obtain 

/ \D(h e - z £ )\ 2 dx < C \h £ - z £ \ 2 dx < d\R £ \ , 

for some positive constants C and C\ which do not depend on e, hence (6.19) holds. 
Applying this argument with 

4 1 

R e — ] — 1,0[ x ] — sin£,sin£[ and R £ = ]—-,—-[ x ]— sine, sine[ 



it follows that (6.17) holds true. 

B-2 




Figure 3. The rectangle R e where we apply a Cacciopoli type estimate 
in order to obtain (6.17). 

From (6.13), (6.14), (6.15), (6.16), and (6.17) we deduce that 

/ \Du £ {a)\ 2 dx < [ \Dv(a)\ 2 dx+ I \Dw £ \ 2 dx + o(l). (6.20) 

JQ e \T(cr) JB- 2 \T((t) JQ £ 



Since 



/ 



Du £ (a)\ 2 dx> I \Dw F \ z dx 



we obtain 

/ \Du £ (a)\ 2 dx < I \Dv(o-)\ 2 dx + o(l) < C (6.21) 

JB_ 2 \r((r) JB- 2 \T(ct) 

uniformly with respect to e. Thus, there exists u*(a) £ H x {B-2 \ r(cr)) such that 

u £ (a) -± u*(a) weakly on H\B- 2 \ r(<r)), (6.22) 

and 

M*(cr) = sin^ on dB_ 2 \ T(a). (6.23) 

As (Du £ (a),D(p) = for every up E H 1 (B_ 2 \ T(a)) with up = on dB_ 2 \ F(a), by 
(6.22) we obtain that (Du*(a),Duj) = 0. By (6.7), this fact, together with (6.23), 
implies that 

u*(a) = v(a) . (6.24) 
In addition, by the lower semicontinuity and by (6.21), we have 

\Dv(a)\ 2 dx < liminf / \Du £ (a)\ 2 dx < I \Dv{a)\ 2 dx. (6.25) 

- 2 \r(a) £ ^ 0+ JB_ 2 \T(cj) JB_ 2 \T{a) 

By (6.22), (6.24), and (6.25), we deduce that (6.8) holds. 



/ 
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